An ultrafast detector for high-power quantum mechanical squeezing 
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Squeezed states of light are a key resource for quantum information, and the detection of squeezing 
is vital in experimental quantum optics. Although squeezed vacuum can be generated at bandwidths 
of 10 THz and more by broadband parametric down conversion, the bandwidth of standard squeezing 
detection methods is limited by the photo-detectors to several GHz at most. We propose sum- 
frequency generation as a physical detector for squeezing, where by using broad phase matching 
techniques, the detection bandwidth can be almost unlimited (up to lOOTHz). We show that by 
measuring the quadrature amplitudes and noise of the sum-frequency light generated by an input of 
squeezed vacuum, the input squeezing level can be accurately deduced, especially with high-power 
input. In addition, the proposed scheme is robust to inefficiency / loss of the sum- frequency photons, 
which reduces the accuracy of the measurement, but does not impede it. 

PACS numbers: 42.50.Dv, 42.50.St, 42.65.Lm, 42.65.Ky 



Quantum mechanical squeezing - the reduction of fluc- 
tuations of one quantum variable at the expense of in- 
creased fluctuations of the conjugate variable - is a ma- 
jor resource in quantum information and quantum mea- 
surement. In optics, squeezed states of light are key to 
methods of phase measurement with precision beyond 
the standard quantum limit (SQL - 1/ \fN , N the total 
number of photons detected), approaching the ultimate 
Heisenberg limit (l/iV) P, Q. Due to the potential for 
dramatic improvement in precision, sub- SQL measure- 
ment methods are appealing for metrology applications, 
such as detection of gravitational waves [3], precision 
spectroscopy 4] and next generation atomic clocks Q . 

A major limitation of standard squeezing detection is 
bandwidth. With homodyne detection, the photo detec- 
tors limit the usable bandwidth to several GHz at most. 
Degenerate parametric down-conversion (PDC) on the 
other hand, can produce ultra broad squeezed vacuum 
states of several tens of THz bandwidth or more, limited 
only by the phase matching bandwidth of the PDC pro- 
cess d, 0, Hi- Thus, a new method is needed to measure 
the quantum properties of ultra-broad squeezed vacuum; 
especially in high power, like those generated by above 
threshold parametric oscillators Uo] . 

Another concern with sub-SQL measurement is the 
sensitivity of squeezed states to photon loss, which cou- 
ples vacuum fluctuations back into the light beam and 
reduces the observed squeezing. Photon loss comes with 
two 'flavors': Loss during propagation and loss due to the 
non-ideal quantum efflciency of photo-detectors. While 
propagation loss can be minimized with optical tech- 
niques (coatings, etc.), the detection loss poses an in- 
herent threshold for the detectability (and usefulness) of 
squeezing with standard homodyne detection. Indeed, 
the best squeezing measurements to date required 
95% — 99% cfBcient detectors, that are hard to come by. 

Sum frequency generation (SFG, also known as para- 
metric up-conversion) was explored previously as a de- 
tector for quantum correlations at very low powers of 



squeezed vacuum input; i.e. individual entangled pho- 
ton pairs [13, [3] ■ III the low power regime, SFG serves 
as an ultrafast two-photon detector, capable of resolving 
simultaneously the tight time-difference and energy-sum 
correlation of the entangled photons. In this letter, we 
examine the quantum properties of SFG at higher powers 
through a full quantum model of SFG light produced by 
squeezed vacuum input of arbitrary power. We introduce 
SFG as a physical detector for high-power squeezed light 
that is both ultrafast and robust to detection inefficiency. 

The detection scheme is outlined in figure [T^ - A local 
oscillator (LO) at frequency 2w, serves as a narrowband 
pump that is down converted to generate squeezed vac- 
uum. The obtained squeezing is detected by up conver- 
sion of the PDC light and homodyning the resulting SFG 
with the original pump LO, while varying the LO phase 
Lp. From the measured quadratures of the SFG light and 
the measured input flux, the squeezing parameters are 
inferred. Interestingly, the scheme is a "symmetric in- 
version" of the standard squeezing detection (outlined 
in figure [lb), where a local oscillator (LO) is first up- 
converted (frequency doubled) and then down-converted 
to generate squeezing, that is later measured by homo- 
dyning with the original LO. In spite of the symmetry, 
the properties of SFG detection are very different, open- 
ing new routs for ultrafast quantum optics. 

Before we review the quantum calculation, it is illu- 
minating to consider a classical model, which although 
simplistic, is very powerful in capturing the main fea- 
tures of its quantum counterpart. Classically, the down- 
converted fleld contains twin frequencies that are corre- 
lated in intensity and anti-correlated in phase {A-^c {uj) — 
A*p(r2p— w)), which in time transforms to a broad- 
band noise fleld with only one quadrature (cdc {t) pure 
real) 0, Loosely speaking, the SFG amplitude 

in time is proportional to the input amplitude squared 
(osFG {t) oc {t) , later we review this assumption thor- 
oughly). Thus, for real input, the SFG output is also real, 
leading to the naive proposition that with squeezed input. 
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FIG. 1: (Color online) (a) - the proposed SFG scheme for 
measurement of squeezing compared to (b) - the standard 
homodyne squeezing detection scheme. 



the SFG output is also squeezed. One must remember 
however that even classically, PDC light is never pure 
real. The PDC process must be seeded by an external 
noise in order to evolve fl^ , and during down-conversion, 
the imaginary quadrature of the seed noise is squeezed, 
but not eliminated. The input PDC field can be written 
as Ojdc 

(t) = q{t) + ip{t), where q,p are the quadrature 
amplitudes. Figure [2^ and [2)3 illustrate the distribution 
of PDC field values O'Dc (^) in the complex plane. Prms is 
reduced by the squeezing ratio M and qrms is increased 
by the same ratio compared to the rms seed noise level hq 
(leaving the product qrms • Prms ~ constant). Quan- 
tum mechanically uq is the vacuum noise nsQL of order 
one photon in amplitude (for convenience we scale our 



units such that 



= 1). The SFG field is then 



OsFG (t) oc (g + ip) = q^ -p^ + 2ipq. 



(1) 



Equation[T]relates the SFG quadratures Q,P to the input 
quadratures q^p, as shown in figures [2j3 and [2ji. The 
imaginary quadrature of the SFG field has zero mean 
with noise of Prms oc qrms ■ Prms « Hg - not Squeezed for 
any M. The real quadrature of the SFG light has non- 
zero average Q ex q'^^^ - pl^^ = - l/AP. For high 
squeezing ratio (M >> 1), where Prms is negligible, the 
PDC flux is a direct measure for 



Q oc $n 



Ir'^ 



Pr 



(2) 



Thus, by measuring the input flux <i>Dc and Prms we 
can infer the input parameters qrms, Prms- With ideal 
squeezing {qrms ■ Prms = 1), the SFG Prms noise will 
be at the SQL level, whereas with non-ideal squeezing 

[irrns ' Prms > 1), thc Prms UOlsC will InCrCaSe. Now 

we can appreciate why SFG detection is robust to loss 
of SFG photons: the squeezed quadrature of the input is 



never directly measured, but rather inferred from the un- 
squeezed noise of the SFG output, so no photo detector 
is required to actually record sub- vacuum fluctuations. 

Let us elaborate further on the properties of SFG with 
broadband input in thick non-linear media. When the 
SFG process is relatively weak (undepleted input approx- 
imation), the SFG field Asfg at frequency Q can be writ- 
ten as a sum over all the contributing frequency pairs in 
the input (down-converted) field A^c (t^) UM- 

Asfg i^) = X2 J dijA^c (w) A^c ~ ij) F (17, Uj) , (3) 

where X2 is the non linear coefficient, assumed to be inde- 
pendent of frequency in the relevant range and F (fi; lo) 
is the phase mismatch function, expressed as 



F (f2, uj) = exp 



.Ak{n,uj)l 



smc 



Ak{n,uj)l 



, (4) 



with / the length of the nonlinear crystal and Ak (fi, oj) 
thc momentum mismatch between the input frequency 
pair {uj,U — uj) and the SFG field at fi. Equation [3] is 
the key to the bandwidth advantage of SFG detection: 
the SFG process conveniently converts the entire input 
bandwidth into one sum-frequency component (the orig- 
inal pump). The broadband squeezing properties can 
therefore be inferred from a narrowband (slow) homo- 
dyne measurement of this one frequency. 

Since we consider SFG as an ultrafast detector, we 
assume broad phase matching conditions, where the mo- 
mentum mismatch can be separated into Ak{Q,Lu) — 
a-{fl-np)+l32-{uj- %/2f + Pi-[uj- np/2f, indicat- 
ing that the output phase matching bandwidth around 
rip is narrow, strongly limited by group velocity mis- 
match (a), whereas the input phase matching bandwidth 
around f2p/2 is much broader, weakly limited by group 
velocity dispersion (/32) and higher orders dispersion {Pn). 
Broad phase matching conditions occur naturally in type- 
I SFG (or PDC) near degeneracy, and especially when 
the center of the down converted spectrum r2p/2 coin- 
cides with the zero dispersion of the nonlinear medium 
{h — 0), where the input phase matching bandwidth can 
span almost an octave, limited only by (3^ [1, [3]. 

The weak (quadratic or higher) dependence of the mo- 
mentum mismatch on the input frequency w, allows to 
neglect it when the input spectrum is properly limited, 
leading to a simple phase mismatch function F that de- 
pends linearly on the output SFG frequency 5^ = fJ — fi„ 



F [Sn) = exp 



.a-l-5n 



a-i-sn 



(5) 



Transformed into time, Eq. [5] represents a rectangular 
phase matching window / {t) = rect [t/to] ; to = al/^i:. 

Consequently, the SFG field amplitude in frequency 
is written as a multiplication Asfg (f^) = X^F {^) x 
J dwAoc (w) ^Dc (f^ — '-^), indicating that the phase 
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FIG. 2: (Color online) Classical simulation of quadrature 
maps for SFG light generated by ideally squeezed input, (a) 
and (b) - input maps for low and high squeezing parame- 
ter M; (c) and (d) - the corresponding SFG output maps 
{Q + iP = {q + ivf)- Note that the real quadrature Q of 
the SFG scales linearly with M^, whereas the imaginary P 
quadrature is independent of M. 



matching function acts as a narrow spectral filter over 
the SFG spectrum. The SFG field in time is therefore 
expressed (using the convolution theorem) as a moving 
window time-average of the square of the input field: 



, (t) = X2 / dt'al,, {t - t') 



(6) 



Note that while the output SFG spectrum is inherently 
filtered by the phase matching bandwidth, it can be even 
further filtered by standard homodyne techniques (mix- 
ing with the original pump). With narrow bandwidth 
homodyning, the averaging time can be extended at will 
beyond t^, up to the total detection time T . 

We now move on to the full quantum calculation of 
the SFG quadratures with squeezed input. We assume 
that the input to the SFG process is a loss free, ideally 
squeezed vacuum state \ipin) — |0, i?) (R being the degree 
of squeezing) In a discrete quantum analog to Eq. 

|6l we write the SFG annihilation operator at time i as a 
sum over input PDC annihilation operators: 



n=T/T 



asFcit) 



4= J2 ^Icit-Jr): 



(7) 



where aoc{t) annihilates a PDC photon in a wave packet 
centered at time t, whose envelope is defined by the 



Fourier limit of the input PDC spectrum. The duration 
of each wave packet is the coherence time r = 1/Adc 
(Adc the input bandwidth), n is the number of wave 
packets summed during the detection time, equal also 
to the number of frequency modes n = T/t = A^c/S 
{S = 1/T is the detection bandwidth). Currently we as- 
sume unit two-photon efficiency for the SFG (x2 — 1 with 
units of field quanta). Later we will remedy this assump- 
tion and the assumption of ideally squeezed input. Fol- 
lowing the standard squeezing calculation method [l7| |. 

a^^. (t) = A{t) cosh R~ (t) sinh R (8) 

where A{t) is the squeezed state annihilator and R is the 
squeezing parameter. We are interested in high squeezing 
{R >> 1), where the squeezing ratio is M — . The 
photon flux of the down-converted squeezed vacuum is 

= ^ (4c(0aDc(t)) = Adc sinh^ R « A^cM^. (9) 

indicating that the input flux is a measure for Af ^ (or R) . 
Using the commutator [A(i), yl"^ (t')] = St^t', we obtain 

a^c (t) = A^{t) cosh^ R + i^' (t) sinh^ R- 
A^t)A{t) + sinh2i?. 



(10) 

Plugging Eq. [10] into Eq. [71 we calculate the SFG photon 
flux [13 $SFG = (asU(i)asFG(i)) /T 

S 



1 



$2 

^DC 



2 

\2 DC 
^DC 



(11) 



Equation[Tl1shows three contributions to the SFG flux - a 
linear term of single pairs (dominating at low squeezing) , 
a coherent term of multiple pairs (2nd, dominating at 
high squeezing), and an incoherent term (3rd, suppressed 
for a broadband input, n = (T/r) >> 1)0, 111, HI- Cal- 
culating for the SFG field the average imaginary quadra- 
ture (PsFG — i (o.sFG{t) — asFG(^))) ^Lud uolsc yields. 



PSFG = 0, (APsfg) = 2. (12) 

Equations [1^ and [5| verify the classical intuition for 
ideally squeezed hght, showing that the SFG imaginary 
quadrature noise is at the SQL level of one (frequency 
doubled) photon in amplitude. We now examine the 
effects of partial squeezing and two-photon inefficiency. 
Loss is generally modeled as a beam splitter (BS) in front 
of an ideal process, as shown in figure[3| In our case, a BS 
at the input of the SFG medium models partial squeezing 
and a BS at the output of the SFG models two-photon in- 
efficiency. Quantum mechanically, these BSs couple vac- 
uum fluctuations into the beams, affecting their noise. 

We begin with the effect of the two-photon efficiency, 
xi is the probability of a single photon-pair to be up- 
converted to an SFG photon. Unity two-photon effi- 
ciency xi — 1 contradicts the undepleted input approx- 
imation. Obviously, the SFG flux cannot exceed the in- 
put flux $sfg < ^Bc, leading to saturation of the SFG 
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FIG. 3: (Color online) Modeling non ideal squeezing / ineffi- 
cient SFG as ideal processes followed by loss beam splitters. 



flux, which for unit two-photon efficiency occurs for any 
R. We note though, that saturation affects mainly the 
real Q quadrature of the SFG, which is non zero and 
large. The imaginary P quadrature is a zero-mean noise 
and is hardly affected. Realistic x| values for common 
non-linear crystals are in the 10~^ — 10~^ range (depend- 
ing on the medium, interaction length, bandwidth, etc.). 
This number can be enhanced several orders of magni- 
tude by resonant cavities [l^ . Accordingly, the transmis- 
sion (reflection) of the output BS in Fig. [3] is ^sfg = X2 
(rgp^_ = 1 — xi)- The flnal output annihilation operator 
can be written as aout = isFcasFc, + irsFcCivac, where avac 
is the annihilator of the vacuum mode entering the BS. It 
is now simple to calculate again the P quadrature noise 
of the final output field, obtaining (APo„t)^ — 2, just like 
the ideal case. Since quadrature noise at the SQL level 
crosses a BS unaffected [3| this result is no surprise. 

The effect of partial squeezing is analyzed in a similar 
manner. When a down-conversion loss BS is introduced 
into the beam before the SFG medium, the input state 
becomes |V'm) = \0,R)i {0)2 ^-nd the annihilation oper- 
ator associated with the input to the SFG medium is 

are the amplitude 

transmission and reflection coefficients of the input loss 
BS. It is now straight forward to plug this input oper- 
ator into Eq. [7] and calculate the imaginary quadrature 
noise again. Using the input photon flux (with the partial 
squeezing loss BS) 



the calculation yields 



{APoutY 



incAnc sinh^ R, 



n 



(13) 



(14) 



Thus, the imaginary quadrature noise increases linearly 
with the input photon loss (r^^,). Equations [T51 and [HI 
allow us to estimate the squeezing parameters M (or R) 
and the input loss r^^ from the measured SFG imagi- 
nary quadrature noise and the known X2- Clearly the 
accuracy of the estimation improves with higher input 
flux (power) and better two-photon efficiency. The two- 
photon (in)cfficiency affects the resolutio'n of the squeez- 
ing estimation, but will not impede the measurement. 



For example, with a 10~^ efficient SFG process, which 
should be accessible with resonant enhancement, the ac- 
curacy of our squeezing estimation will be reduced by a 
factor of 10* compared to the ideal SFG detector, but 
will still yield a valuable result. For comparison, trying 
to detect squeezing with 10~* efficient photo detectors 
using standard homodyne techniques is futile. 

In summary, a new concept has been proposed for 
measuring quantum mechanical squeezing, using SFG as 
an ultrafast, physical correlation detector. This method 
opens a new avenue to measuring the quantum properties 
of broadband, high power squeezed states, as produced 
by parametric oscillators. The physical correlation detec- 
tion by the two-photon SFG process makes the proposed 
method robust to detection inefficiency, as opposed to 
standard homodyne techniques. The new method can be 
applied to applications of precision metrology, such as 
ultra-precise interferometry for gravitational waves de- 
tection, as will be discussed in future publications 
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